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Abstract 

Thermal properties in the metal/graphene (Gr) systems are analyzed by using an atomistic 
phonon transport model based on Landauer formalism and first-principles calculations. The spe- 
cific structures under investigation include chemisorbed Ni(lll)/Gr, physisorbed Cu(lll)/Gr and 
Au(lll)/Gr, as well as Pd(lll)/Gr with intermediate characteristics. Calculated results illustrate 
a strong dependence of thermal transfer on the details of interfacial microstructures. In particular, 
it is shown that the chemisorbed case provides a generally smaller interfacial thermal resistance 
than the physisorbed due to the stronger bonding. However, our calculation also indicates that the 
weakly chemisorbed interface of Pd/Gr may be an exception, with the largest thermal resistance 
among the considered. Further examination of the electrostatic potential and interatomic force 
constants reveal that the mixed bonding force between the Pd and C atoms results in incomplete 
hybridization of Pd and graphene orbital states at the junction, leading effectively to two phonon 
interfaces and a larger than expected thermal resistance. Comparison with available experimental 
data shows good agreement. The result clearly suggests the feasibility of phonon engineering for 
thermal property optimization at the interface. 

PACS numbers: 63.22.-m,65.80.-g,73.40.Ns,66.70.-f 
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I. INTRODUCTION 



As the critical dimensions of modern electronic devices approach the nanoscale, the power 
density of integrated circuits has soared drastically, forcing the issue of thermal manage- 
ment to the forefront. Faced with the challenge, graphene has recently emerged as one of 
the key candidates for the next-generation low-power electronics for its superb properties. 
Fundamental understanding of thermal transport in the graphene based structures is cru- 
cial from the perspective of both low-dimensional physics and practical applications of this 
emerging material system. Particularly, the metal contact with graphene is of interest as 
it provides not only an essential part of any active device but also a primary path of heat 
dissipation. Thermal conduction across the heterogeneous metal/graphene (Gr) interface 
is characterized by the interfacial resistance often known as the Kapitza resistance.!^ The 
heat current is mainly carried by phonons, while the electronic contribution is consider- 
ably smaller. It was found recently that thermal energy transfer via direct elect ron-phonon 
coupling astride the heterointerface is insignificant (compared to the phonon-phonon inter- 
actions) in metal/dielectric structures including those employing graphene.^ 

As in the case of most layered structures, the properties of the metal/graphene system 
are also expected to depend heavily on the bonding chemistry and detailed structures at the 
interface. It was recently demonstrated that the bonding between graphene and metal atoms 
can be divided into two categories— chemisorption and physisorption.^ While chemisorption 
opens a band gap in graphene due to hybridization between the graphene p^- and metal 
d-orbitals, the Dirac-cone feature of graphene is preserved at the physisorbed interfaces. 
Subsequent investigations illustrated that physisorption is observed for Ag, Al, Cu, Cd, Ir, 
Pt, and Au, whereas the Ni, Co, Ru, Pd, and Ti interfaces belong to chemisorption.!^!^ 
However, their impact on phonon/thermal transport has received much less attention with 
only a limited number of studies available in the literature. A molecular dynamics calculation 
was conducted for heat transfer between allotropic carbon nanotubes and Cu substrate.^ On 
the side of experiment, measurement of thermal conductance was reported for the Cu/Gr, 
Al/Gr, Ti/Gr and Au/Gr cases with the obtained values ranging 0.8 — 2.5 x 10~^ Km^/W^ 

In this paper, we attempt to provide a detailed theoretical account of interfacial thermal 
resistance in the metal/graphene system. The sample structures are chosen to refiect the 
range of typical interfaces from chemisoprtion to physisorption for a comprehensive analy- 
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sis. Since atomistic details of phonon dynamics are crucial for the accurate outcome, the 
adopted theoretical approach utilizes a first-principles analysis based on density functional 
theory (DFT) to calculate the interatomic force constant (IFC) matrix.flSHn] xhen the de- 
sired phonon transmission function and the thermal current are determined via Green's 
function technique^^^^ and the Landauer formalism.^ The obtained results are compared 
with experimental and other theoretical data available in the literature. 

The rest of the paper is organized as follows. First, the structures of metal/graphene 
systems under investigation are discussed briefiy with the focus on the interfacial bonding 
environment. Then, a summary description of the adopted thermal transport model is 
provided including the details of numerical implementation. Finally, discussion/analysis of 
the calculated interfacial phonon transmission and Kapitza resistance is presented. 

II. METAL/GRAPHENE INTERFACES 

The heterogeneous system of interest is graphene on the (111) surface of a metal (Ni, Pd, 
Cu, and Au). If the lattice constant of graphene is fixed at the optimized value 2.46 A, less 
than 5% lattice mismatch is introduced when these metals of face-centered cubic symmetry 
are made commensurate with the graphene lattice. A 1x1 unit cell is formed for Ni/Gr and 
Cu/Gr, while a 2x2 construction is necessary for Pd/Gr and Au/Gr to accommodate larger 
sizes of metal atoms. The top view and stacking orders of the four diflFerent material combi- 
nations are shown schematically in Fig. [T| representing the ideal atomic arrangement at the 
metal/graphene contacts. While useful in illustrating the geometric construction, these ar- 
rangements are insufficient to account for details of the interfacial microstructures. Instead, 
the realistic structures can be obtained through geometry optimization using first-principles 
calculation (see, for example. Fig. [2]). As indicated earlier, the relaxed metal/graphene inter- 
faces are generally divided into two categories, chemisorption and physisorption, depending 
on the bonding energies, interfacial separations, and orbital hybridizations. In accord with 
earlier studies,!^ our calculation clearly illustrates that Ni and graphene bond strongly at the 
interface through hybridization between Ni rf-orbitals and C p^-orbitals. Strong coupling 
also results in a relatively small interfacial separation (2.02 A). On the other hand, Cu and 
Au are physisorbed on graphene and form a weak van der Waals bonding with a larger 
interlayer distance (2.89 A and 3.31 A, respectively).^^ For the Pd(lll)/Gr structure, it 
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is observed that the interaction between Pd and C atoms may be smaUer than the strong 
chemical bonding, leading effectively to the mixed character at the interface (i.e., of both 
chemisorption and physisorption).^ This unique combination may set the Pd/Gr system 
apart from the other metal/graphene interfaces as elaborated further in Sec. IV. 



III. THEORETICAL MODEL 



A. Thermal transport in the atomistic Green's function formulism 

Phonon transport is originated from the dynamics of the lattice or lattice vibrations.!^ 
Since the investigation concerns the interfacial properties, transport in the immediate region 
astride the interface can be treated ballistic and a quantum mechanical treatment in the 
Landauer framework!^ is adequate to capture the essential features. Accordingly, we consider 
a three parted system where the central interface region (i.e., the region of interest) is 
connected to the thermal reservoirs on the left and the right with two semi-infinite leads 
(labeled L and R), often known as the lead-conductor-lead configuration.'^^^ Accounting for 
only the phonon contribution, the thermal current density is then given as: 



Jin,Tn) = —nuTA'^)[n{n,u)-n{Tn,uj)], 



(1) 



where n(TL,R, uj) is the equilibrium phonon distribution in the left or right thermal reservoir 
(at temperature Tl,r) and Tph{^) denotes the phonon transmission function. In this expres- 
sion, Tph{^) is directly related to the lattice dynamics of the given structure containing all 
of the relevant details. Its computation can be achieved in a manner analogous to that of 
electron transmission coefficients across nanoscale heterostructures by using such approaches 
as the atomistic Green's function met hod. '^^'^^'^ 

The first step of calculating the thermal current [i.e., Tph{^)] is to obtain the IPCs defined 

as: 

K,,,^A^) = ^^^ for ^ 7^ J, (2) 

where E is the total energy of the system, R is a Bravais lattice vector, and Ui^a is the dis- 
placement of the i*'^ atom (of the unit cell) in the a direction with respect to the equilibrium 
position. The i = j cases can be determined in terms of those in Eq. ^ via the acoustic 
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sum rule.'^ The IFCs in turn define the harmonic matrix 



and the dynamical matrix 



(3) 



(4) 



where Mi is the mass of the i*^ atom. In the considered lead-conductor-lead system, this 
harmonic matrix can be written formally o^^^ 



K = 



( i^LC ^ 



kl^ kc k^c 
kL k 



(5) 



RC / 



and, finally, the dynamical equation as 



det 



= 0. 



(6) 



Here, Kp is the matrix in each of the three regions (p = L,C,R), while Kpc represents the 
phonon coupling at the interfaces between the central region and the two leads {p = L,R). 
This expression in Eq. ([g]) is clearly analogous to the governing equation of the electronic 
system, with uj'^ ^ e and K ^ H [see, for example, Eq. (1) of Ref. [15]. Accordingly, calcu- 
lation of phonon transmission can directly follow the Green function techniques developed 
initially for electronic transport as mentioned earlier. 

Utilizing the identified parallelism, the phonon transmission function is written straight- 
forwardly as: 

nu) = Tr {V^G^X^G}.) , (7) 

where G\f are the retarded (r) and advanced (a) Green's functions of the central region 
and rL,R correspond to the coupling with the left and right leads, respectively. Further, the 
Green's function and the coupling functions can be obtained frompH 

-1 



C 



{u + iQY - Kc 



where the self-energy terms 



= ^Ic^pkpc 



(8) 
(9) 

(10) 



are evaluated with the help of the Green's function in the corresponding leads {p = 
L,R). The transfer matrix method offers a very efficient approach to calculate in the 
semi-infinite lead by treating it as a stack of principal layers with only nearest-layer inter- 
actions .^^'^^'211 xhe advanced terms (with superscript a) are given as Hermitian conjugates. 
Additional details of underlying theoretical formulation can be found in Ref. [21] 

Once the phonon transmission function T{uj) is obtained, thermal and thermoelectric 
properties can be evaluated. A particularly relevant quantity is the thermal resistance Rth 
of the system that is given as the inverse of the thermal conductance: 

duj dn{T^(jjy ^ 



f 

Jo 



While the model described above is for the case of two leads (or reservoirs), extension to a 
multi-terminal system is trivial as it has been demonstrated in electronic transport.'^ 



B. Numerical implementation 

The calculations are performed in the DFT framework, as it is implemented in the 
QUANTUM-ESPRESSO package,'^ with ultrasoft pseudopotentials in the local density 
approximation.^^ A minimum of 35 Ry is used for the energy cut-off in the plane wave 
expansion along with the charge truncation of 350 Ry. In addition, the Methfessel-Paxton 
first-order spreading is employed with the smearing width of 0.01 eV. The momentum space 
is sampled on a Monkhorst-Pack mesh in the first Brillouin zone. In the construction of 
interfacial structure, two layers of graphene and three to five layers of metal (depending on 
the material) are considered in the calculation. The unit cell is set to the lattice constant 
of graphene and the geometry optimization is performed to find the energy minimum struc- 
ture. The Hellman-Feynman force is taken into account in each iterative solution. Figure [2] 
shows the resulting interface structures of the considered material combinations. They serve 
as the central region in the previously mentioned lead-conductor-lead configuration. Two 
leads consisting of respective bulk materials (i.e., bulk metal and graphene/graphite) con- 
nect seamlessly to the interface region and are modeled separately. No appreciable mismatch 
(i.e., resistance) exists between the leads and the conductor. Simulation results show that 
the k-space grids of 6x6x2 and 6x6x6 are sufficient for the interface region and the bulk 
leads, respectively, with good convergence in the relevant characteristics. 
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With the optimized reahstic interface structures, the next step is to evaluate vibrational 
properties including IFCs. While analytical force constant models have been used widely 
in such materials as carbon or silicon based systems,'^^''^ they are not directly applicable 
in the current investigation due to the lack of information on the required bonding details 
between carbon and metal atoms. Instead, we calculate lattice dynamics, also within the 
DFT formalism for accuracy, by utilizing a perturbative treatment known as density func- 
tional perturbation theory (DFPT).'^The IFCs are then obtained by Fourier analyzing the 
set of dynamical matrices generated from the first-principles calculation under the harmonic 
approximation [see Eq. With the IFCs of bulk leads, we can build K^^ {p = L,R) 

that enter Eq. ([H]); While the on-site matrix Kq and coupling matrices K^c and Krc can 
be constructed from the IFCs of the central part. Then the phonon transmission can be 
evaluated using Eq. ([7| with the Green's function calculated from the transfer matrix tech- 
nique. Finally, the thermal resistance of the structure and the Kapitza resistance at the 



metal/graphene junction can be evaluated by Eq. (11). Since this quantity of interest is 
defined in the near-equilibrium condition (i.e., an infinitesimally small temperature gradient 



across the structure; see Eq. (11)], an equilibrium treatment is adequate with no need for 
an iterative solution. 



IV. RESULTS AND DISCUSSION 

Of the material combinations under investigation, graphene on the Cu (111) surface is 
examined first due, partly, to its widespread use. As indicated by the obtained phonon 
transmission function in Fig. [3](a), only the low lying acoustic branches (below 100 cm~-^) 
play the dominant role in phonon transport at the interface. The impact of optical branches 
is orders of magnitude smaller. The resulting thermal resistance is plotted in Fig. IsF 



Since the total resistance of the structure [i.e., Eq. (11)] contains the contribution from 
the leads as well, the intrinsic thermal resistance at the junction is deduced by subtracting 
this portion in a manner analogous to electrical transport J ^^ ^ ^^ * At T=300 K, the estimated 
interfacial thermal resistance of the Cu/Gr structure is 1.18 x 10~^ Km^/W. The interfacial 
resistance exhibits the 1/T dependence in the low temperature region (50—150 K), while 
staying almost invariant between 150 K and 450 K. 

When the calculation is extended to chemisorbed Ni and physisorbed Au, the respective 



7 



results appear to be similar in many aspects. However, one interesting point to note is that 
the Ni/Gr interface shows the phonon transmission coefficient whose frequency dependence 
is much broader with fewer resonant features (see Fig. [3]). This is substantially different from 
those of the physisorbed metal/graphene interfaces (both Au and Cu). The mixed nature of 
phonon dynamics at the chemisorbed interface (thus, a smaller mismatch) is thought to be 
the main origin of enhanced phonon transmission and eventually a smaller interfacial resis- 
tance. The estimated value for Ni/Gr is about 3.9 x 10~^ Km^/W at 300K, whereas it is more 
than four times larger for Au/Gr (1.7 x 10~^ Km^/W) as indicated in Fig. [5] The obtained 
result for physisorbed Au/Gr shows good agreement with the measurement data available 
in the literature (^ 2 — 3 x 10~^ Km^/W).'^In contrast, the conventional diffuse mismatch 
model shows large discrepancies. Even with the added sophistication such as anisotropy in 
graphitic materials and multiple heat transfer mechanisms, it substantially overestimates 
the thermal resistance for Au/Gr (^ 6.8 x 10~^ Km^/W).*^ A similar calculation yields 
- 4 X 10"^ KmVW at the Cu/Gr interface.EH 

Comparison between different graphene/metal systems thus far clearly indicates that the 
chemisorbed interface is generally more favorable than the physisorbed in term of thermal 
transport. The presence of strong bonding and the smaller interlayer separation (e.g., 2.02 
A of Ni/Gr vs. 3.31 A of Au/Gr) all support this conclusion that is also in accord with a 
recent experimental study.^ One potential exception may be graphene on the Pd (111) sur- 
face. Since their bonding characteristics supposedly show both chemisorbed and physisorbed 
nature as mentioned earlier, it is reasonable to anticipate that the interfacial thermal re- 
sistance would fall in between as well. However, the calculation suggests that Pd may not 
follow the trend and actually have the largest resistance of those considered (a value of 
3.35 X 10~^ Km^/W at room temperature; see Fig.[6]). Most notably, phonon transmission is 
greatly suppressed between approx. 40 cm~^ and 100 cm~^ even when compared to Cu/Gr 
and Au/Gr. This result is counterintuitive, particularly when the interlayer distance (which 
tends to indicate the interaction strength) behaves as expected; namely, between the values 
of chemisorbed and physisorbed structures as indicated in Fig. [2j 

The identified peculiarity is examined further by analyzing detailed microstructures of 
the corresponding metal/graphene interfaces. The electrostatic potential isosurface plotted 
in Fig. [7| clearly illustrates that strong hybridization between the and d orbitals at the 
chemisorbed Ni/Gr interface "glues" Ni and graphene together, making them essentially 
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one unit. As such, propagating phonons transmit through the Ni/Gr interface with relative 
ease. On the other hand, the weak interaction between Au and C atoms in the Au/Gr 
system forms a barrier at the interface, which strongly scatters phonon transmission. For 
Pd/Gr, it is revealed that the mixed bonding force between Pd and C atoms is indeed 
smaller than the strong chemical bonding (namely, chemisorption) as evident from the partly 
detached bond in the boxed region. However, the interaction still alters the orbital states of 
first layer graphene, making them sufficiently distinct from those of second layer graphene. 
Consequently, phonons potentially face two interfaces for transmission instead of one. A 
force constant analysis between atomic layers can clarify the latter point with numerical 
certainty. 

Figure [8] illustrates the interlayer force constants deduced from the DFPT calculation. 
The height of each bar symbolizes the interaction strength between two neighboring layers. 
For example, the first bar on the left denotes the interaction between layers 1 and 2; the next 
bars are for layers 2 and 3, and so on. In all three plots, the metallic layers are up to layer 
5 and graphene starts from layer 6; accordingly, the physical interface of two heterogeneous 
materials is located between layers 5 and 6. On the other hand, the real interface or barrier 
which impedes phonon transport is characterized by the abrupt change of force constant. 
In Ni/Gr, it is illustrating to note that the force constant between the Ni and graphene 
layers right at the interface (the bar between layers 5 and 6) shows only a slight difference 
with those between Ni layers on the left (i.e., 0.103 a.u. vs. 0.118 a.u., where a.u. stands for 
atomic Rydberg units). Instead, transmitting phonons experience the major barrier at the 
interface between layers 6 and 7, where the force constant changes the most drastically. Due 
to strong hybridization discussed earlier, the first layer of graphene is absorbed by Ni atoms, 
leaving it practically decoupled from the second graphene layer. Nevertheless this is a Gr/Gr 
interface and it is reasonable to expect a relatively smaller resistance. For physisorbed Au, 
the characteristics are as expected. Namely, we get a clear separation of gold bonding and 
graphene bonding at the interface between layers 5 and 6 (0.028 a.u. vs. 0.004 a.u.). Since 
gold atoms are much heavier than carbon, acoustic phonon frequencies experience a large 
mismatch — hence, a large thermal resistance. When it comes to Pd/Gr, the incomplete 
mixing of the first graphene layer indeed results in two-step changes in the force constant 
between layers 5 and 6 as well as between 6 and 7 (i.e., from 0.04 a.u. to 0.0125 a.u. 
then to 8.31 x 10""^ a.u.), leading effectively to two phonon interfaces as mentioned earlier. 
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and a larger than expected interfacial thermal resistance. One cautionary point is that the 
calculation outcome could experience modifications if the graphene film is just one monolayer 
thick. Then, the thermal resistance values can be substantially smaller than the presented, 
particularly when significant mixing is involved. Our current theoretical formalism is not 
equipped to address isolated systems, for which the ideal leads are difficult to construct. 

Additional insight into interfacial phonon transport may be gained by comparing the re- 
sults of metal/graphene structures with those involving dielectric substrates. As summarized 
in Table|lp2'^a couple of likenesses can be readily noted. That is, the calculated interfacial 
thermal resistance of Ni/Gr is close to the corresponding value of BN/Gr (3.9 x 10~^ vs. 
5.4 X 10-^ KmVW), while Pd/Gr and SiC/Gr are very much alike (3.35 x 10"^ vs. 3.61 x 10"^ 
Km^/W). Since it is the Gr/Gr interface in Ni/Gr that determines the thermal resistance, 
the chemisorbed case can be understood roughly analogous to the epitaxial Gr/BN structure, 
also consisting of two materials of similar two-dimensional (2D) crystal type. In the case of 
Pd/Gr vs. SiC/Gr, both harbor effectively more than one interface/barrier for phonon trans- 
mission. As it is well known, the studied structure of epitaxial graphene grown on the SiC 
(0001) surface contains an additional carbon buffer layer at the interface that does not have 
the characteristic sp^ bonding. Hence, the observed large interfacial thermal resistance in 
SiC/Gr is consistent with the discussion on Pd/Gr given earlier (i.e., two phonon barriers in 
series) . Finally, the physisorbed structures with a weak van der Waals bonding between two 
dissimilar materials (e.g., 3D metal vs. 2D graphene) may correspond to exfoliated graphene 
placed on a non-2D crystal substrate such as Si02. While comparable theoretical estimate 
based on a first principles calculation is not available for Si02/Gr, the recent experimen- 
tal data^^ of 0.56 — 1.2 x 10~^ Km^/W show good match with the values of Cu/Gr and 
Au/Gr obtained in the current investigation. These analyses highlight how significantly the 
atomic bonding details can infiuence the interfacial thermal properties, leading potentially 
to phonon engineering for active heat management .^^^ It also illustrates the inadequacy 
of various theoretical treatments including the molecular dynamics approachP^ that cannot 
properly account for the required level of physics a priori. 
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V. SUMMARY 



Thermal transport in the metal/graphene heterostructures is investigated by using a first 
principles method and the Green's function approach within the Landauer formalism. The 
obtained interfacial thermal resistances are 3.9 x 10~^ Km^/W, 1.18 x 10~^ Km^/W, and 
1.70 X 10~^ Km^/W at room temperature for the Ni/Gr, Cu/Gr and Au/Gr structures, 
respectively, indicating generally more effective thermal transfer at the chemisorbed surface 
owing to the smaller interlayer separation and stronger bonding with graphene. However, 
calculations also illustrate that a weakly chemisorbed case such as Pd/Gr could actually 
lead to an interface even more resistive than that encountered at the physisorbed surface. 
Detailed examination of electrostatic potential and force constants identifies the formation 
of an intermediate layer (a consequence of incomplete mixing) and the resulting multiple 
phonon interfaces/barriers as the potential origin of the observed deviation. Comparison 
with the corresponding calculations in the graphene/substrate systems reveals strong corre- 
lation between seeming differences in material combination, further emphasizing the role of 
atomic-level ab initio analysis. The obtained theoretical results show good agreement with 
experimental data available in the literature. 
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TABLE I: Thermal properties at the relevant graphene/metal and graphene/substrate interfaces. 
The entries for Si02/Gr are from experiments, whereas the other are from the first-principles 
calculation. 





Interface 
environment 


Interfacial 
separation 


Thermal resistance 
(10-10 W/W) 


Ni/Gr 


Chemisorption 


2.02 A 


39 


Cu/Gr 


Physisorption 


2.89 A 


118 


Au/Gr 


Physisorption 


3.31 A 


170 


Pd/Gr 


Mixed 


2.43 A 


335 


BN/Gr-^ 


Flat 


3.43 A 


54 


SiC/Gr« 


Rough (buffer layer) 


3.89 A 


361 


SiOs/Gr'' 


Rough 


4.2 A 


56-120 



Ref. M 
Ref. El 
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FIG. 1: (Color online) Schematic representation of (a) stacking order and (b,c) top view of graphene 
absorbed on the metal (111) surface, (b) corresponds to Ni/Gr and Cu/Gr, while (c) is for Pd/Gr 
and Au/Gr. Di specifies the interfacial separation between graphene and the metal. 
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FIG. 2: (Color online) Side view of the investigated metal/graphene systems; (a) Ni(lll)/Gr, 
Pd(lll)/Gr, (c) Cu(lll)/Gr, and (d) Au(lll)/Gr. Two upper layers represent graphene. 
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FIG. 3: (Color online) (a) Phonon transmission function vs. frequency and (b) interfacial thermal 
resistance vs. temperature for Cu/Gr. The vertical dashed line marks the resistance at room 
temperature (300 K). 
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FIG. 4: (Color online) Phonon transmission function vs. frequency for Ni/Gr and Au/Gr. 
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FIG. 5: (Color online) Interfacial thermal resistances vs. temperature for Ni/Gr and Au/Gr. The 
vertical dashed lines mark the resistances at room temperature (300 K). 



19 



140 



360 




0.0 0.5 1.0 50 150250350450 
Transmission Function Temperature (K) 



FIG. 6: (Color online) (a) Phonon transmission function vs. frequency and (b) interfacial thermal 
resistance vs. temperature for Pd/Gr. The vertical dashed line marks the resistance at room 
temperature (300 K). 
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FIG. 7: (Color online) Electrostatic potential isosurface of the metal/graphene system for (a) 
Ni/Gr, (b) Pd/Gr, and (c) Au/Gr. Top three layers are for the metal, while bottom two correspond 
to graphene. In (b), partly detached bonding in the boxed region indicates incomplete mixing at 
the Pd/Gr interface. The large (colored) balls denote the metallic atoms. 
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FIG. 8: (Color online) Interlayer force constants for Ni/Gr, Au/Gr and Pd/Gr. The height of each 
bar represents the interaction strength between two layers, where a.u. stands for atomic Rydberg 
units. In all three plots, the metallic layers are up to layer 5 (i.e., 1—5) and graphene starts from 
layer 6 (i.e., 6—10). The dashed lines indicate the physical interface between metal and graphene, 
whereas the dash-dotted lines (as well as the block arrows) symbolize the barriers that transmitting 
phonons experience. 
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